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Although carbon nanotubes consist of honeycomb carbon, they have never been fabricated 
from graphene directly. Here we show, by quantum molecular-dynamics simulations and classical 
continuum-elasticity modeling, that graphene nanoribbons can, indeed, be transformed into carbon 
nanotubes by means of twisting. The chiralities of the tubes thus fabricated can be not only pre- 
dicted but also externally controlled. This twisting route is a new opportunity for nano-fabrication, 
and is easily generalizable to ribbons made of other planar nanomaterials. 
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I. INTRODUCTION 

Carbon nanotubes (CNTs) are today grown by atomic 
self-organization, resulting in variable tube diameters 
and structures 12 , whereas in textbooks nanotubes are 
portrayed as being made by "rolling- up graphene". Re- 
cent experiments have shown how nanotubes can be un- 
zipped into graphene nanoribbons (GNRs)P^ The in- 
verse route of making nanotubes of graphene directly, 
however, has remained elusive. 

Here we show that graphene nanoribbons can indeed 
be transformed into carbon nanotubes by twisting. Our 
quantum molecular-dynamics simulations show that tube 
formation is preceded by buckling that serves as an in- 
termediate stage between a flat ribbon and a pristine 
tube. Both buckling and tube formation can be explained 
by classical continuum elasticity, apart from effects due 
to atomic discreteness. Since graphene nanoribbons can 
be made with controlled edge morphology 5 , even with 
atomic precision and sub-nanometer width 6 , this twist- 
ing route to carbon nanotube fabrication provides a qual- 
itatively new opportunity for precision control in nan- 
otube fabrication, while also enabling encapsulation of 
molecules, designed chirality, and easy generalization to 
ribbons made of other planar nanomaterials. 



II. THE TUBE FORMATION PROCESS: 
SHOWCASE 

To demonstrate the mechanism of tube formation, we 
first present a simulation of the twisting of a 24 A wide, 
infinitely long graphene nanoribbon with unpassivated 
zigzag edges, shown in Fig. la. We simulated the fixed- 
length ribbon with quantum molecular-dynamics at room 
temperature, while increasing the twist at a constant 
rate; increasing the twist continuously without end effects 
was enabled by periodic boundary conditions adapted to 
chiral symmetry (see Appendix [B] for simulation details). 
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The amount of twist is characterized by a dimension- 
less parameter rW with r being the twist angle per unit 
length and W the width; hence rW = for a flat ribbon 
and rW = 1 for a ribbon that has been twisted a full 




Twist rW 

FIG. 1: (Color online) Quantum molecular-dynamics simula- 
tion of GNR transformation into a CNT. (a) The simulation 
of a 12-ZGNR with snapshots at different twist parameters 
rW (for GNR notations, see Appendix |A|. The red atoms 
highlight one row of atoms across the ribbon, and their cross- 
section projections are shown below; the cross denotes the 
twist axis. The rightmost tube is the final, fully-optimized 
(9, 3) CNT. (b) The tubularity parameter Q that shows the 
stages of flat ribbon, buckled ribbon, and CNT. Fluctuations 
in Q arise from non-zero temperature, (c) Potential energy 
per unit length as a function of twist. The final curve seg- 
ments in panels b and c (the green ones) represent the subse- 
quent full structural optimization of the resulting CNT. 
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turn within a length of 2ttW. 

At the initial phases of the simulation, the cross sec- 
tion of the ribbon remained flat until a critical value 
rW = TbW = 0.83 when the ribbon buckled into a 
twisted groove with a U-shaped cross section, see Fig. la. 
In order to characterize the tubular geometry quantita- 
tively, we introduced a geometrical parameter, the tubu- 
larity parameter which measures how much the ribbon 
edges have approached each other in relation to ribbon 
width: it is zero for flat ribbons and one for tubes [see 
Eq.(|AlJ]. 

Upon further twisting, ribbon's tubularity was fluctu- 
ating, but increased on the average. Increased twisting 
brought therefore the opposite edges closer together, and 
at a second critical value, rW = r t W = 1.62, they be- 
gan to interact chemically: they were joined by a sudden 
formation of bonds, and the buckled ribbon was rapidly 
transformed into a tube. This point was identified as a 
sudden increase in the tubularity parameter Q and as an 
even more sudden decrease in the potential energy due 
to the formation of a and tt bonds. At this stage, the 
resulting tube had a slightly non-round cross section. So 
as to anneal any remnant residual strains, we stopped 
the molecular-dynamics simulation and carried out a full 
structural optimization. The simulation then resulted in 
a pristine (9,3) CNT (see Supplemental Video l). We 
repeated such simulations for a number of zigzag and 
armchair nanoribbons of varying width, and invariably 
obtained pristine CNTs. 

The essence of the tube formation process outlined 
above can be captured by a simple demonstration — just 
twist the ends of a strap of your backpack and watch the 
result (Fig. 2). In particular, the CNT formation process 
is not an artifact of the imposed boundary conditions; 
simulations with finite tubes and thousands of atoms pro- 



duce the same results, as discussed in Section VI On the 



actual nano-scale, twisting experiments should be fea- 
sible using the established paddle-type setups in which 
voltages applied to electrodes can be used to control the 
orientation of the somewhat asymmetrically positioned 
paddle electrostatically (the inset in Fig. 2a). This setup 




FIG. 2: (Color online) Ribbon buckling and tube formation 
can be demonstrated with a strap of a backpack. Inset: Pro- 
posal for an experimental realization, adapting the successful 
scheme of Refs [8] and |9] 



has already been successful in twisting CNTsP^ 

These atomistic simulations showed the feasibility of 
tube formation by twisting, but raised several questions 
that include characterization of buckling and tube forma- 
tion as a function of ribbon width, the required critical 
torque, the resulting CNT chiralities and their possible 
control, electronic structure modifications, finite-size ef- 
fects, effects of imperfections, and, finally, a simple ex- 
planation of the mechanisms involved at a classical con- 
tinuum level. We attempt to address these questions in 
the following sections. 



III. TUBE FORMATION IS GOVERNED BY 
CONTINUUM ELASTICITY 

We addressed the last question by modeling GNRs as 
thin elastic sheets with an in-plane modulus fc, bending 
modulus K, Poisson ratio v = 0.3, and width W. The 
twisted shape was found by numerically minimizing the 
elastic deformation energy that had an in-plane stretch- 
ing component and an out-of-plane bending component 
(see Appendix[C|. Our analysis based on continuum elas- 
ticity theory indicates that the key elements in tube for- 
mation are the following: Buckling results from trans- 
verse stress in the twisted ribbon and the shape of the 
buckled cross section from competition between stretch- 
ing and bending. 

The condition of fixed length in the analysis distin- 
guishes it from previous studies in which the tensile force 
was fixed, and ribbon lengt h cou ld thus vary leading to 
smooth longitudinal buckling 10 11 , helical developable ge- 
ometry 12 , or triangular stress-focusing patterns^. Our 
analysis assumed that the cross section of the ribbon 
was free to bend and warp, which is not possible close 
to its ends if they are clamped. The ribbon should thus 
be sufficiently long for our assumptions to be valid. The 
length 27t/t£, over which the ribbon makes a full turn at 
the point of tube formation, can be taken as a reasonable 
minimum-length criterion. 

Dimensional analysis indicates that the twisted geom- 
etry depends only on the dimensionless elastic thick- 
ness, h = \/K/ (kW 2 ), and twist, rW. For graphene 
K = 1.6 eV and k = 25 eV/A 2 such that h w 0.01 
in the above simulation with W ~ 24 Figure 3b 

shows the changing geometry of an elastic ribbon un- 
der twisting for h — 0.01 — the similarity with Fig. [iji 
is evident. Quantitative comparison in Fig. (3}d of the 
tubularity parameters of quantum and classical analyses 
makes the similarity of these two approaches even more 
evident, although there are some differences near Q w 1. 
This is because of the difference between the zigzag and 
armchair ribbons of similar width, which have somewhat 
different tube-formation mechanisms due to the role of 
edge morphology and edge stress^M that do not nor- 
mally appear in continuum models (discussed in Section 
fV| ). We note here that the maximum strain induced by 
twisting at the ribbon edge is « 6% for h = 0.01, decreas- 
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FIG. 3: (Color online) Tube formation by elasticity theory, (a) Visualizations of twisted elastic ribbons for h = 0.01 (width 
w 24 A in graphene) and rW — 0, 0.44, 0.88, 1.32, and 1.65. Cross-section projections are also shown, (b) Tubularity parameter 
Q by elasticity theory (h — 0.01) and from simulation results for a zigzag (12-ZGNR, h ~ 0.0104) and an armchair ribbon 
(22-AGNR, h ~ 0.0097). (c) Buckling (r&) and tube-formation points (r t ) as a function of elastic thickness, which result from 
elasticity theory and from simulations for different ZGNRs. For the rightmost ribbon, 10-ZGNR, the simulation was repeated 
10 times, and the vertical bars denote the range where t& and r t were observed to fluctuate, (d) Buckling and tube-formation 
points that result from elasticity theory (h — 0.012) and from simulation results for 10-ZGNR (h ~ 0.0117), as a function of 
pre-strain. 



ing as 1.3 A/W for increasing W, so graphene ribbons 
survive twisting without tearing. 



The buckling point, t&, and the tube formation point, 
rt, vary with the scaled elastic thickness, h, as shown 
in Fig. [3]3. These numerically obtained curves are well 

approximated by the expressions r^W « 6.6 yh and 

r t W « Y7.2yh that follow from a simple scaling anal- 
ysis (see Appendix [C]). The same analysis also yields 
the torque associated with the tube-formation point as 

M t « l.SK/y/h. In the presence of an externally applied 
pre-strain, 70, which naturally can be released after tube 
formation, the critical values for 75 and r t decrease as 
shown in Fig.[3]i. This is in agreement with the trend ob- 
served in molecular-dynamics simulations: Control over 
70 thus provides a way to fine-tune the required twist and 
has implications for the resulting CNT chiralities, which 
is the question we discuss next. 



IV. CNT CHIRALITIES CAN BE PREDICTED 

Our simulations suggest that, by knowing the width 
and chirality of the GNR, the chirality of the CNT can 
be predicted with unexpected reliability. The CNT chi- 
ral angle, depends on a shift, Az, that measures the 
(relative) axial displacement of ribbon edges at the tube- 
formation point. In the continuum picture the honey- 
comb geometry thus implies the relation 

(j) + = it/2 - arctan(Az/W) (1) 

between the GNR (0) and CNT (0) chiral angles^ 3 that 
vary between 0° and 30° (see Appendix [A|. With Az/W 
given by the elasticity theory and with tube circumfer- 
ence deduced from the ribbon width (Fig. [9]), Eq.Q offers 
a recipe for a continuum prediction of the CNT chirality. 

In the atomistic picture, however, the shift Az has 
to be compatible with the atomic discreteness. Fig- 
ure [4^i shows how Az in Af-ZGNRs picked values close 
to the continuum predictions, either an integer or a half- 
odd-integer times the edge periodicity a zz = 2.46 A. 
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FIG. 4: (Color online) Predicting CNT chiralities. (a) Axial shift Az of the opposite edges in the tube as given by the elasticity 
theory as a function of ribbon width (lines) for three values of pre-strain (denoted by percentages). Circles are the simulated 
values for ZGNRs. (b) The CNT chiral angles, 6, that correspond to the shifts in panel a. Full lines are the predictions for 
given by Eq.(T| with = 0° and with Az/W given by the elasticity theory, (c) Shift Az for 10-, 12- and 14-ZGNRs 
(symbols from left to right) and the corresponding elasticity-theory prediction (full curve) with zero pre-strain. Each of these 
ribbons was simulated 10 times, and the intensity of the color inside the circles is proportional to the number of times the 
room-temperature simulation gave the given Az. (d) The CNT chiral angles, augmented by the chiral indices, corresponding 
to the Az of c, with the same definitions for the symbols. The line is the elasticity-theory prediction by Eq.|l]). 



Other values of would have cost additional shear- 
deformation energy. Similarly, Fig. [4}} shows that the chi- 
ral angles of the corresponding CNTs' pick allowed values 
in the proximity of the continuum-limit curves given by 
Eq.0. 

The elasticity result for Az can, in fact, be used for a 
still better estimate of the CNT chiral angle — also for an 
estimate of the CNT chiral indices (n, m) themselves. For 
a given N and ribbon type, only certain indices (n, m) 
are allowed, and the honeycomb lattice suggests the ex- 
pressions 



then given by (n, m) via the expression 



(n, m) 



\N/2^S h ,N/2- 
(N/2 + 5\N/2- 
(N/2-5 h ,25 h ), 

{{N/2-8\2V), 



- 5 h ), odd-TV ZGNR 
5% even-TV ZGNR 
odd-TV AGNR 
even- TV AGNR, 



(2) 



where 8 % is the nearest integer and 5 h is the nearest half- 
odd-integer to the ratio Az/a[ zz / ac ], owing to the sym- 
metry of the opposite edges as shown in Fig. [8] (TV/2 can 
be a half-odd-integer) , and where Az is predicted by the 
elasticity theory. Prediction for the CNT chiral angle is 



0(n, m) 



arccos 



2n + m 



2Vr, 



(3) 



Yet Eq.Q is only a prediction — the tube formation 
process contains stochastic aspects. For some ribbons 
the continuum- limit value for Az happened to be halfway 
between two allowed values, making prediction based 
on room-temperature simulations inevitably less precise. 
For example, Fig. and [4]i show how identical room 
temperature simulations result in different shifts Az and 
in correspondingly different CNT chiralities. While the 
obtained shifts indeed favored the nearest integer multi- 
ple of a zz close to the elasticity-theory prediction, they 
fluctuated so that the final result could not be predicted 
with certainty. Still for more than 80 % of our simula- 
tions Eq.Q predicted the CNT chirality correctly. 

Finally, since increasing the axial pre-strain 70 de- 
creases rt, and thus Az, as shown in Fig. |4^, we can use 
70 to control the CNT chirality as illustrated in Fig. ^p. 
In particular, given the discreteness of the allowed Az, 
even a coarse experimental control over 70 can be used 
to fine-tune the chirality. 



5 




Elastic thickness h{xlO 2 ) wjdth w width W (A) 



FIG. 5: (Color online) Tube formation for armchair graphene nanoribbons. (a) Buckling (rfc) and tube-formation points (rt) as 
a function of elastic thickness given by the elasticity theory (full lines; Eq.(IT]) with (j> — 30°) and by simulations (squares) for 
AGNRs of varying width (corresponds to Fig. ^ for ZGNRs). The inset: AGNR prior to tube formation, illustrating the radial 
bulging of the edges. Pre-strain is zero, b) Shift Az as given by the elasticity theory (full lines) and by atomic simulations 
(squares) as a function of ribbon width for three values of pre-strain (corresponds to Fig. [4^l for AGNRs) (c) The CNT chiral 
angles corresponding to the shifts in b. 



V. ARMCHAIR RIBBONS: EDGE STRESS AND 
FORCED-JOINING EFFECTS 

Figure (5^i shows the buckling and tube-formation 
points for AGNRs, showing less apparent agreement be- 
tween simulations and continuum-limit theory than for 
ZGNRs (compare with Fig. [3^). These differences have 
two reasons: compressive edge stress and 'forced-joining 
effect' that arises from the comparatively large unit-cell 
length at the ribbon edge. 

The compressive stress at the edge of an unpassivated 
AGNRs is^ 1.45 eV/A, some ~ 3.5 times greater than in 
ZGNRsP In narrow ribbons this stress can cause spon- 
taneous twisting because of elongation of the edge with 
respect to the ribbon axis. 18 Because the stress makes 
edges to prefer small strain, larger twists were required 
for buckling (consistent shifts between the blue curve and 
the blue symbols in Fig. [5^l) . However, we omitted the 
edge stress from the elasticity theory on three grounds. 
First, it affected mainly the buckling threshold. Second, 
the edge stress depends on edge type and passivation — 
with hydrogen passivation the stress vanishes. Third, the 
edge stress can be imitated simply by having a decreased 
pre-strain. 

Apart from buckling, the tube formation itself was 
dominated by a kind of forced-joining effect caused by 
edge morphology. Namely, because the edge periodicity 
a ac = 4.26 A in AGNRs is almost twice the periodicity 
a zz = 2.46 A in ZGNRs, larger deviations from the con- 
tinuum prediction for Az were required — certain 'forcing' 
was needed to initiate the tube formation. Joining was 
easier for wide ribbons, where shifts in steps of a ac re- 
quired less shear, and for ZGNRs, where the required 
steps a zz were smaller. 

Furthermore, prior to tube formation the edge stress 
turned our to make the edges of buckled ribbons to bulge 
radially outwards (the inset of Fig. Then, to initi- 
ate the tube formation, the buckled ribbon with bulged 



edges often required more forced twisting so as to attain 
the allowed Az. In some simulations bulging was further 
enhanced by twisting, making tube formation to require 
exceptionally large twists (tube formation was hindered 
by an energy barrier). In these simulations a longer sim- 
ulation time or higher temperature might have initiated 
the tube formation earlier. However, when the contin- 
uum Az happened to be such that dangling bonds from 
the opposite edges met directly, no forcing was needed 
and simulation yielded r t in agreement with the elastic- 
ity theory prediction (compare Figs[5|i and^). 

Note that while pre-strain drives ZGNRs towards arm- 
chair CNTs (0 = 30°; see Fig.[4]o), it drives AGNRs to- 
wards zigzag CNTs (0 = 0°; see Fig. [5J3). These ten- 
dencies can be understood by noticing that increased 
pre-strain always decreases Az; with precisely zero Az 
armchair ribbons would become zigzag tubes and zigzag 
ribbons armchair tubes. 

In addition to the purely geometrical effects associ- 
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FIG. 6: (Color online) Energy gaps for selected N- AGNRs 
from quasi-static simulations (optimized atomic coordinates 
for given rW), representing three different N families. The 
narrow panel shows the gaps in the resulting CNTs with 
torque and axial stress removed. 



FIG. 7: (Color online) Finite ribbon simulations. Figure shows snapshots from the simulation of L — 38 nm long 12-ZGNR with 
clamped ends and with increasing end-to-end twist angle (the first six snapshots; structures are rotated to give the best view 
angles). The last snapshot is the (partially unzipped) (9,3) CNT after removing the clamps and performing a constraint-free 
simulation at 300 K. 



ated with the GNR-CNT transition, we also investigated 
twisting-induced modifications in the electronic proper- 
ties of the system. Prior to buckling, the energy gaps 
in iV-AGNRs were observed to form three "families" ac- 
cording to q = mod(7V, 3), as reported earliei^Ej, anc [ 
as illustrated in Fig. [6j Buckling turned out to cause 
only slight rehybridization of the carbon atoms, consis- 
tent with earlier studies on CNTs. 22 In contrast to the 
smooth electronic modification caused by twisting, at the 
point of tube formation the gap instantaneously jumped 
to that of CNT — which may have been, prior to relax- 
ation, still influenced by remnant residual torques. 

VI. NOT AN ARTIFACT OF 
PERIODIC BOUNDARY CONDITIONS 

The validity of the revised periodic boundary condition 
(RPBC) approach was confirmed by finite-ribbon simula- 
tions. Here we present exemplary results of a simulation 
for L = 38 nm long 12-ZGNR with 3720 atoms. Rib- 
bon's one end was clamped, and the other end was kept 
at a fixed distance (without pre-strain) while twisted con- 
tinuously at the rate dO/dt — 0.0091 degrees/fs with a 
2 fs time step; the rest of the atoms were treated with a 
Langevin thermostat set to 300 K. In this simulation we 
used the REBO interatomic potential from the LAMMPS 
packagePSES 

Supplementary Video 2 shows an animation of this 
simulation, with selected snapshots in Fig. [7| The tube- 
formation processes in RPBC and in finite ribbons were 
the same — and both resulted in the same pristine (9, 3) 
CNT. Sure enough, some finite-size effects did arise. 
Both buckling and tube formation initiated in a narrow 
central region, and the zipping-up propagated towards 
the ribbon ends when the applied twist increased. Com- 



plex distortions were suppressed by the experimentally 
feasible fixed- length constraint. After tube formation we 
released the end constraints and observed that the tubes 
(that were partially unzipped near the ends) remained 
thermodynamically stable at 300 K, and even at 1200 K. 

The spreading of the end-to-end twist angle AO across 
the ribbon was somewhat uneven, and twisting took ef- 
fectively place within a length smaller than L. Therefore, 
while initiation of the buckling at r^W = {A6\ ) /L)W ~ 
0.83 agreed with what happened for RPBC, the initiation 
of tube formation at r t W = (A0 t /L)W « 0.91 occured 
earlier than in RPBC. For increasing length such finite- 
size effects vanished and the results converged towards 
RPBC results; for large L buckling and tube- formation 
points also became sharper. We performed this simula- 
tion four times for different twist rates, obtaining invari- 
ably the same results. 

Such trends in the finite-size effects were confirmed by 
simulations of shorter ribbons. For instance, dimensional 

analysis helped to find the scaling L/W > 0.7^/w/A as 
the critical length-to- width ratio above which the picture 
of the tube-formation process remains valid. Indeed, for 
W = 25 A, inferring a critical ratio of ~ 4, tube was 
formed in an expected manner for L/W = 7.5, but not 
for L/W = 2.6. More systematic investigations of the 
finite-size effects are underway. 



VII. CONCLUDING DISCUSSION 

When the twisted ribbons have atomically smooth 
edges, which is experimentally feasible and even 
preferred^, the formed CNTs are exp ected to become 
essentially pristine by energy arguments! 26 * 27 ! In practice, 
however, we cannot exclude the formation of defects ei- 
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ther. If tube formation is initiated at different locations 
with different CNT chiralities 28 , the zipping- up of the 
tube may give rise to scattered point defects. Moreover, 
edge roughness, irregular edge chirality, and edge passiva- 
tion can lead to CNTs with vacancies, impurity atoms, or 
dangling bonds, arranged as chiral line defects. Although 
we cannot entirely exclude the appearance of phenomena 
related to other finite-size effects^, lattice fatigue^, or 
complex defect formation^, preliminary results indicate 
that the central concepts of tube formation prevail. Fur- 
thermore, when GNRs are hydrogen-passivated, as they 
often are, tube formation must be preceded by dehydro- 
genation and formation of H2 gas. Since this reaction has 
only a weak thermodynamic driving force^, presumable 
energy barriers for formation of H2 suggest a slow reac- 
tion, and catalytic dehydrogenation may be required to 
aid the tube formation.^ 

To conclude, our study opens up new opportunities in 
nanomaterial manipulation not limited to carbon-based 
ribbons alone. Indeed, using a combination of varying ge- 
ometry that ribbons afford with their separation of scales, 
one might envisage using inhomogeneous width, chemical 
modification including passivation, adsorbed molecules, 
and clusters, to construct structures with new function- 
alities. Examples include tubes with bulges or partial 
tears^, nanoscrolls, multiwalled nanotubes with a spi- 
ral cross section^, all of which can also be manipulated 
using external forces so as to enable molecular encapsu- 
lation and release in a variety of applications. 
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Appendix A: About ribbon geometries 

We simulated tube formation for 7V-ZGNRs (0 = 0°, 
N = 10 - 19) and AGNRs (0 = 30°, N = 16 - 34), 
their geometries are illustrated in Fig. [8j All simulated 
ZGNRs and AGNRs resulted in pristine CNTs with well- 
defined chiral indices (n, m) (corresponding to CNTs 
uniquely defined by the vector C = na + mb, the cir- 
cumferential vector expressed in terms of the honeycomb 
unit-cell vectors a and b). For computational feasibility 
ribbons were unpassivated; hydrogen passivation would 
have required catalyst particles or a prohibitively long 
simulation time. 

The tubularity parameter was defined as 
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FIG. 8: (Color online) ZGNRs and AGNRs illustrating the N 
alternation in the opposite-edge profiles; alternation shows up 
e.g. in Eq.(2]). The unit-cell length is a zz « 2.46 A in ZGNRs 
and a ac ~ 4.26 A in AGNRs for zero pre-strain. Ribbon 
relaxation leads to slightly longer ribbons, more so for narrow 
ribbons. Width Wo is measured from atomic positions [Wo ~ 
(2.13 • N - 1.42) A for ZGNRs and W « 1.23 • (N - 1) A for 
AGNRs]. The angle (f) is the angle between ribbon axis and 
the adjacent zigzag direction. 



where (ibond = 1-42 A is the carbon-carbon bond length 
and d is the distance between any two opposite-edge 
atoms that form a bond in the final tube. The distance is 
at maximum (d = dfl a t) f° r a flat ribbon and at minimum 
(d = ^bond) fc> r a tube. Hence Q = for a flat ribbon and 
Q = 1 for a tube. The threshold for buckling was defined 
as ^buckled > 0.1. In the continuum limit, because there 
are no bonds, Eq.(Al) was used with (ibond 



0. 



n = 



^flat 



^flat — ^bond 



(Al) 



The width of an atomistic ribbon is a question of def- 
inition, and although the width Wo (Fig. [8]) would be 
an easy concept, the direct comparison of atomic and 
continuum widths is inherently ambiguous. We chose to 
define the ribbon width W of the atomistic ribbon as the 
circumference of the resulting CNT (Fig. [9J. The differ- 
ence between W and Wo, which mainly originates from 
curvature and bond formation (merging of opposite edges 
creates 'new surface area'), plays a bigger role in narrow 
ribbons (when dbond is a notable fraction of Wo)- 
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FIG. 9: (Color online) The ratio of CNT circumference to 
Wo- Fitting gave for the tube circumference ttD = W ~ Wo + 
2.5 A (Poisson effect caused by the pre- strain was removed). 
The CNT circumference was obtained more accurately from 
the chiral indices, W = v^bondV^ 2 + m 2 + mn with (n, m) 
predicted by Eq.|2|. 



Appendix B: Molecular-dynamics simulations 

We used spin-unpolarized density-functional tight- 
binding^ESl anc [ rev i se d p eriodic bo undary conditions 
adapted to chiral symmetry! 18 * 21 * 37 ^ In the twisting sim- 
ulations, with minimal cell in the axial direction and zero 
strain corresponding to a relaxed flat ribbon, we used the 
Langevin thermostat at 300 K with 1 fs time step, and 
a stepped twist rate of A(rW)/At = 0.2 - 0.3 ns -1 . 
(Buckling and tube formation was possible due to the 
absence of symmetry constraints with respect to axial 
symmetry, unlike in Refs fTTJI and I2T1) The rate has only 
a minor effect on the results because of the abruptness 
of the buckling and tube-formation events. Molecular- 
dynamics simulations were performed for ZGNRs using 
20 ^-points and for AGNRs using 10 ^-points (while cal- 
culating energy gaps using 100 ^-points) with respect to 
the chiral symmetry operation. 



Appendix C: Analysis of twisting of a ribbon based 
on elasticity theory 



We consider twisting of a thin ribbon with a fixed 
length and translational symmetry such that each cross 
section of the ribbon has the same shape. A cross sec- 
tion is free to warp in the direction of the twist axis. We 
denote by x = (xi^x 2 ) the material coordinates of the 
ribbon, where x\ and x 2 are the coordinates in the trans- 
verse and longitudinal directions, respectively. Position 
in space of a material point x is given by r(x) = (x, y, z), 
and the twist axis was chosen to coincide with the z-axis. 



Symmetry of the problem implies that 

dr 

dx 2 

and 



(-ry,rx,l + 70) 



d 2 r 

dx\ 



(~T 2 X, -T 2 ?/,0), 



(CI) 



(C2) 



where r is the twist per unit length and 70 is a longitu- 
dinal external strain. In-plane deformation of the sheet 
is described by the strain tensor 



1 / dr dr 

lij = 2\dx~ i 'dx~ j ~ ij 



(C3) 



and out-of-plane deformation by the curvature tensor 



c d 2 r 

13 dxidxj 



(C4) 



where n is the surface normal and z,j = 1,2. Deforma- 
tion energy per unit length is given b}^ 



U 



k r w 

2(1 + 1/) Jo 



Tr 7 2 



1 - 2v x 



(Tr 7 ) 2 Uci (C5) 



K 



/ { (Tr C) 2 + (1 - v) [Tr C 2 - (Tr C) 2 ] } dx x 
Jo 



Here k is the in-plane modulus, K the bending modulus, 
v the Poisson ratio and W the ribbon width. The shape 
of the ribbon was found by numerically minimizing U. 
To this end we discretized the cross section x 2 = const 
into N = 100 points and replaced the derivatives by finite 
differences. The discretized energy with v = 0.3 was min- 
imized by a damped molecular-dynamics method. Both 
75 and r t could be determined by increasing r in small 
steps. 

Useful insight into the buckling and tube formation can 
be obtained by a simple scaling analysis. During twist 
the initially straight longitudinal 'fibers' (narrow strips 
across the ribbon) are deformed into helices and strained 
by ^T 2 i? 2 , where R is the distance from the twist axis. 
This generates a compressive stress T ~ kr A W A in the 
transverse direction of a flat ribbon as the fibers tend 
towards the axis to minimize the longitudinal stretch- 
ing energy. The buckling threshold for a compressed 
plate is given 

bjJSI T - K/W 2 , from which a critical 

twist, r b W « c b Vh, is obtained. Here h = ^/K/(kW 2 ) 
is the elastic thickness. For r > 75 the buckled shape 
of the cross section is determined by competition be- 
tween stretching and bending. The stretching and bend- 
ing energies per unit length can be estimated such that 



5 « § (\t 2 R 2 ) W and B 
cross section curved with a radius R. By minimizing 



KW 
2R 2 ' 



respectively, for a 



S + B for R 



w 



2?r , we find that the twist r t W w c t Vh 
is required to bring the two edges of the ribbon to- 
gether so as to form a tube. For the torque we find that 
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FIG. 10: (Color online) (a) Torque at the tube- formation point as determined by numerical energy minimization (full line) is 
compared with the scaling estimate Mt/K ~ 1.8/ \/h (dashed line), (b) Tube-formation point, rt, and buckling point, t&, are 
compared with the scaling estimates rt ~ 17. 2y/h and r b w 6.6V% respectively. 



M t 



d(S+B) 



cmK/ v ft. The constants C5 = 6.6, 



q = 17.2, and cm = 1.8 were found by fitting results of 
numerical energy minimization, see Fig. [To] The simple 
scaling expressions work well, although deviations appear 
at high ft when the cross-section warps are large. 

During twisting nearly all of the shear strain vanishes, 
which leads to warping of the cross section, i.e., to rela- 
tive displacement of the two edges along the twist axis. 
Integrated warp, or the shift Az, can be approximated 
by 



Az 



W/2 
W/2 



rr ( • — ) dx 



where r is the distance from the twist axis and is a unit 
vector perpendicular to the radial direction. By assuming 
a circular cross section with radius R = ^- , we find 



Az 



T t W 2 

2tt 



(C7) 



(C6) for the shift at the tube-formation point. 



' email: pekka.koskinen@iki.fi 

1 S. Iijima, Nature 354, 56 (1991). 

2 R. Saito, G. Dresselhaus, and M. S. Dresselhaus, Physical 
Properties of Carbon Nanotubes (Imperial College Press, 
London, 1998), 1st ed. 

3 D. V. Kosynkin, A. L. Higginbotham, A. Sinitskii, J. R. 
Lomeda, A. Dimiev, B. K. Prince, and J. M. Tour, Nature 
458, 872 (2009). 

4 L. Jiao, L. Zhang, X. Wang, G. Diankov, and H. Dai, Na- 
ture 458, 877 (2009). 

5 X. Jia, M. Hofmann, V. Meunier, B. G. Sumpter, 
J. Campos-Delgado, J. M. Romo-Herrera, H. Son, Y.-P. 
Hsieh, A. Reina, J. Kong, et al., Science 323, 1701 (2009). 

6 J. Cai, P. Ruffieux, R. Jafaar, M. Bieri, T. Braun, 
S. Blankenburg, M. Muoth, A. P. Seitsonen, M. Saleh, 
X. Feng, et al., Nature 466, 470 (2010). 

7 N. Hamada, S. Sawada, and A. Oshiyama, Phys. Rev. Lett. 
68, 1579 (1992). 

8 A. M. Fennimore, T. D. Yuzvinsky, W.-Q. Han, M. S. 
Fuhrer, J. Cumings, and A. Zettl, Nature 424, 408 (2003). 

9 J. C. Meyer, M. Paillet, and S. Roth, Science 309, 1539 



(2005) . 

10 A. E. Green, Proc. R. Soc. London 154, 430 (1936). 

11 A. E. Green, Proc. R. Soc. London 161, 197 (1937). 

12 E. H. Mansfield, The Bending and Stretching of Plates 
(Cambridge University Press, Cambridge, 1989), 2nd ed. 

13 A. P. Korte, E. L. Starostin, and G. H. M. van der Heijden, 
Proceedings of the Royal Society of Longon. Series A 47, 
285 (2010). 

14 P. Koskinen and O. O. Kit, Phys. Rev. B 81, 235420 
(2010). 

15 H. Liang and L. Mahadevan, Proc. Nat. Ac. Sci. 106, 22049 

(2009) . 

16 B. Huang, M. Liu, N. Su, J. Wu, W. Duan, B. Gu, and 
F. Liu, Phys. Rev. Lett. 102, 166404 (2009). 

17 V. Barone, O. Hod, and G. E. Scuseria, Nano Lett. 6, 2748 

(2006) . 

18 P. Koskinen and O. O. Kit, Phys. Rev. Lett. 105, 106401 

(2010) . 

19 D. Gunlycke, J. Li, J. W. Mintmire, and C. T. White, 
Nano Lett. 10, 3638 (2010). 

20 D.-B. Zhang and T. Dumitrica, Small 7, 1023 (2011). 



10 



21 P. Koskinen, Appl. Phys. Lett. 99, 013105 (2011). 

22 V. N. Popov, New J. Phys. 6, 17 (2004). 

23 



S. Plimpton, J. Comp. Phys. 117, 1 (1995). 
http : //lammps . sandia . gov. 



24 

25 D. W. Brenner, O. A. Shenderova, J. A. Harrison, S. J. 
Stuart, B. Ni, and S. B. Sinnott, J. Phys.: Condens. Matter 
14, 783 (2002). 

26 S. Malola, H. Hakkinen, and P. Koskinen, Phys. Rev. B 
81, 165447 (2010). 

27 P. Koskinen, S. Malola, and H. Hakkinen, Phys. Rev. Lett. 
101, 115502 (2008). 

28 S. Cranford and M. J. Buehler, Modelling Simul. Mater. 
Sci. Eng. 19, 054003 (2011). 

29 K. V. Bets and B. I. Yakobson, Nano Res 2, 161 (2009). 

30 M. B. Nardelli, B. I. Yakobson, and J. Bernholc, Phys. 
Rev. Lett. 81, 4656 (1998). 

31 Y. Li, J. Phys. D: Appl. Phys. 43, 495405 (2010). 



32 N. Shah, Y. Wang, D. Panjala, and G. P. Huffman, Energy 
& Fuels 18, 727 (2004). 

33 C. Tang, W. Guo, and C. Chen, Phys. Rev. B 83, 075410 
(2011). 

34 L. M. Viculis, J. J. Mack, and R. B. Kaner, Science 299, 
1361 (2003). 

35 D. Porezag, T. Frauenheim, T. Kohler, G. Seifert, and 
R. Kaschner, Phys. Rev. B 51, 12947 (1995). 

36 P. Koskinen and V. Makinen, Comput. Mater. Sci. 47, 237 
(2009). 

37 T. Dumitrica and R. D. James, J. Mech. Phys. Solids 55, 
2206 (2007). 

38 O. O. Kit, L. Pastewka, and P. Koskinen, Phys. Rev. B 
84, 155431 (2011). 

39 L. D. Landau and E. M. Lifshitz, Theory of elasticity 
(Pergamon, New York, 1986), 3rd ed. 



